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Purpose: Primary congenital glaucoma (PCG) is an autosomal recessive form of glaucoma that manifests within the first 
year of life and if left untreated, leads to irreversible blindness. Cytochrome P450 1B1 (CYP1B1) is the major gene known 
to be associated with PCG. The role of the CYP1B1 gene in disease pathogenesis and the relatively low detection rate of 
CYP1B1 mutations in some populations, especially Asians, remain unexplained. We hypothesized that altered gene dosage 
of CYP1B1 or anterior segmental dysgenesis causative genes may be involved in the pathogenesis of PCG. 
Methods: We performed whole genome exon-focused array comparative genome hybridization (aCGH) to identify copy 
number variation (CNV) in 20 Korean PCG patients and their parents. 

Results: We identified 12 patients with at least one rare gene-containing copy number variation each, corresponding to 
25 CNVs (5 deletions and 20 duplications) at frequencies of 5-30% in PCG patients and 0% in controls. The 25 CNVs 
were not located at known chromosomal loci for PCG, namely GLC3A, which harbors CYP1B1 (2p21), GLC3B (lp36.2- 
p36.1), or GLC3C (14q23), and did not include any target genes associated with PCG or anterior segmental dysgenesis. 
Conclusions: Further genetic studies with larger cohorts of patients are necessary to validate our results and to elucidate 
other genetic mechanisms underlying PCG, because the identified CNVs might be PCG-specific pathogenic variants and 
may explain the disease pathogenesis of PCG. 



Primary congenital glaucoma (PCG), an autosomal 
recessive disorder of the eye that usually manifests within the 
first year of life, is caused by developmental defects in the 
trabecular meshwork and anterior chamber angle [1]. These 
developmental defects cause obstruction of aqueous outflow, 
leading to increased intraocular pressure and optic nerve 
damage. If left untreated, PCG leads to irreversible blindness 
[2]. 

Three loci for PCG (gene symbol GLC3), namely 
GLC3A (chromosome 2p2-p21), GLC3B (chromosome 
Ip36.2-p36.1), and GLC3C (chromosome 14q23), have been 
mapped by linkage analysis [3-5]. Cytochrome p450 1B1 
(CYP1B1- OMIM 601771), which is associated with GLC3A, 
is known to be the major gene associated with PCG. The 
coding and promoter regions of CYP1B1 have been screened 
extensively in PCG patients worldwide. However, the 
pathway by which CYP1B1 affects development of the 
anterior chamber of the eye remains unknown. In addition, the 
prevalence of CYP1B1 mutations varies in different 
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populations, ranging from -10% in Mexico to 100% in 
consanguineous Slovakian Gypsy patients. The incidence of 
CYP1B1 gene involvement in Korean, Chinese, and Japanese 
PCG patients is relatively low (-20%), which means that it is 
not possible to explain the genetic pathogenesis of PCG in 
more than half of patients who do not have any CYP1B1 
mutations or only one mutant CYP1B1 allele [6-10]. To try to 
explain these findings, we noted that anterior segmental 
dysgenesis (ASD) disorders such as Axenfeld-Rieger 
syndrome, Peters' anomaly, and iris hypoplasia, are 
associated with raised intraocular pressure and an increased 
incidence of glaucoma. Furthermore, the eye is exquisitely 
sensitive to both reduced and increased doses of key 
developmental genes such as paired box 6 (PAX6) and 
forkhead box CI (FOXC1), as demonstrated in ASD [11-15]. 

In this context, we expanded the disease spectrum of PCG 
to ASD and hypothesized that identification of underlying 
genomic imbalances could lead to elucidation of other genetic 
mechanisms underlying PCG. Therefore, our aim in this study 
was to explore alternative genetic mechanisms related to 
disease pathogenesis in PCG. We performed whole genome 
exon-focused array CGH (aCGH) to investigate the most 
susceptible copy number variations (CNVs) in 20 PCG 
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Table 1. Clinical data including CYP1B1 and MYOC sequencing results of subjects with primary congenital glaucoma. 



Individual 
number 

1 

2 
3 
4 
5 
6 



9 

10 
11 
12 
13 
14 
15 
16 
17 



18 
19 
20 



Gender 

M 
M 
M 

F 

F 
M 

F 

F 

F 
M 

F 

F 

F 

F 

F 
M 
M 

M 
M 
F 



Age of onset 
(months) 

10 

3 

76 
21 
<1 
29 

<1 
<1 
117 

5 

3 

1 

35 
10 
5 
6 
4 

2 
10 
1 



Affected 
eye(s) 

U 
B 
B 
B 
B 
U 

U 

u 

B 

u 
u 
u 
u 

B 

u 
u 

B 

B 

U 
B 



IOP at diagnosis 
(OD;OS) 

22 
22/16 
30/24 
34/17 
30/32 

34 

23 
22.4 
11/44 
33/15 

28 
18.5 

30 
52/52 

28 
22.5 
37/35 

26/32 
18.5 
21/30 



Family CYP1B1 mutation MYOC mutation 
history 



Sibling 



L228S 
(heterozygote) 



A330F 
(heterozygote) 



Abbreviation: U, unilateral; B, bilateral. 



patients and their parents. To the best of our knowledge, this 
is the first gene dosage analysis of PCG patients using aCGH. 

METHODS 

Patients and clinical evaluation: Twenty Korean PCG 
patients and their unaffected parents were recruited during a 
multi-institutional collaborative study from September 2008 
to February 2010 and included in this study. Clinical data 
including CYP1B1 and myocillin (MYOC) sequencing results 
of subjects with PCG are shown in Table 1. Criteria for PCG 
diagnosis were evaluated when examination was possible and 
included IOP >21 mmHg in at least one eye; megalocornea; 
corneal edema/clouding/opacity; and glaucomatous optic 
nerve head damage. Corroborating features included 
symptoms of epiphora and photophobia. Patients with other 
ocular or systemic anomalies were excluded. Of the 20 PCG 
patients, only two patients had a known heterozygous 
mutation and a novel mutation in CYP1B1 and MYOC, 
respectively. 

Array CGH analysis: To identify highly susceptible CNVs in 
PCG compared with control individuals, we performed aCGH 
on 20 PCG patients and their parents along with 99 healthy 
individuals. Genomic DNA was extracted from peripheral 
blood leukocytes using a Wizard Genomic DNA Purification 
kit according to the manufacturer's instructions (Promega, 
Madison, WI). A total of 159 DNA samples were labeled and 
co-hybridized to determine DNA copy number changes in 
deletion/duplication using the NimbleGen Human CGH 
3x720K Whole Genome Exon-Focused Array (Roche 
Diagnostics, Mannheim, Germany). Arrays were washed and 



then scanned using a NimbleGen MS 200 Microarray Scanner 
with 2 jum scanning resolution. Raw copy number data were 
normalized using Nexus software 5.0 (Nexus BioDiscovery, 
El Segundo, CA). All probe coordinates were mapped to the 
National Center for Biotechnology Information (NCBI) 
human genome assembly Build 36 (UCSC hgl8). The 
normalized data were then processed with quality controls 
using Nexus Software and the manufacturer's recommended 
default settings. To remove wave artifacts, we checked the 
number of aligned probes with signal quality and applied 
lowess correction to the log2 ratios. The same CNV types were 
merged with adjacent CNV calls using the criteria of <10 
probes and <50 kb apart. The data were loaded into the Fast 
Adaptive States Segmentation Technique (FAS ST) 
segmentation algorithm with a significance threshold of 
1 .Ox 10~ 5 . To minimize the number of false positive CNV calls 
without compromising the sensitivity of detection of true 
CNVs, we applied log2 ratios with thresholds of 0.3 in gain 
signals and -0.5 in loss signals. We visually inspected each 
sample, and excluded CNVs with <5 probes or CNVs <500 
bp in length. 

Identification of copy number variations: The CNV filtering 
strategy that we used is summarized in Figure 1 . To describe 
the candidate CNVs, we used the Database of Genomic 
Variants (DGV; The Centre for Applied Genomics, Toronto, 
Canada) to determine whether the CNVs were novel or 
known. We excluded known CNVs in DGV and those 
identified in controls, and determined candidate genes using 
the identified CNVs. 
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20 PCG patients and their parents (60 samples) 
99 healthy controls without glaucoma phototype 



Whole genome exon-focused array CGH (NimbleGen, Roche, Germany) 



Nexus software (all quality controls were conducted according to the manufacturer's 
recommendations such as Log^ratios, coverage, depth, length 500 bp, and _5 probes 



Deletions 
(PCG patients: 13 CNVs) 



Duplications 
(PCG patients: 142 CNVs) 



Filter 1: Exclusion of known CNY regions based on the Database of Genomic Variants 



Filter 2: Exclusion of CNY regions identified in the 99 healthy controls 



Identification of 25 candidate CNVs (28 genes) in PCG patients 
(PCG patients: 5 deletions and 20 duplications) 



Type of CNV 


Analysis of inheritance mode 

25 candidate CNVs detected in 20 PCG patients matched against their 
parents' CNVs, 


Loss of 1 copy 


One copy loss from one parent = Detected in 5 CNVs 


Loss of 2 copies 


Not detected 


Gain of 1 copy 


One cop)' gain front one parent (a) = Detected in 9 CNVs 


De novo variation (b) - Detected in 6 CNVs 


Mixed cases with (a) and (b) = Detected in 5 CNVs 


Gain of 2 copies 


Not detected 





Filter 3 or 4 

Filter 3: Selection of CNVs located in known chromosomal loci for PCG: GLC3A harboring 
CYP1B1 (2p21), GLC3B (lp36,2-p36,lX and GLC3C (14q23) 



Filter 4: Selection of CNVs compatible with target genes associated with PCG or ASD or 
their highly predicted functional partner genes 



The most susceptible CNVs (genes) in PCG patients 



Figure 1 . Summary of the copy number variation filtering strategy used in this study. 
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Table 2. Table of predicted functional partners of the CYP1B1, MYOC, LTBP2, PAX6, PITX2, PITX3, FOXC1, FOXE3, 

EYA1, LMX1B, and MAF genes by STRING 8.3. 



CYP1B1 

COMT 

GSTM1 

GSTP1 

UGT1A1 

EPHX1 

GSTM3 

UGTB7 

UGT1A6 

GSTA1 

UGT1A9 

FOXC1 

PITX2 

DLL4 

HEY2 

ALX4 

MSX2 

PAX6 

FGF2 

FLNA 

CYP1B1 

CXCR4 



MYOC 

OPTN 

CYP1B1 

WDR36 

SNCG 

FN1 

SERPINF1 

HRAS 

TMTC1 

TIMP1 

OLFM3 

FOXE3 

GPR161 

SOX2 

SMAD4 

PITX2 

FOXF2 

FOXH1 

SMAD2 

HDAC9 

POU5F1 

RSRFR2 



LTBP2 

LTB 

TGFBI 

ZNF185 

ELN 

LOXL1 

MFAP2 

COMMD9 

ESYT3 

ACYP1 



EYA1 

SIX1 

SIX2 

DA CHI 

PAX2 

GDNF 

SIX5 

PAX1 

TLX1 

MYOG 

PAX3 



PAX6 

SIX3 

MITF 

SOX2 

NEUROG2 

SOX3 

IP013 

PKNPX1 

TRIM 11 

SHH 

GCG 

LMX1B 

NPHS2 

WNT7A 

NPS 

COL4A4 

PRODH2 

GPM6A 

COK4A3 

FOXA2 

LDB1 

LRP6 



PITX2 

CTNNB1 

SUCLG1 

LEF1 

FOXC1 

CCND2 

GIPC1 

RGS1 

RGS2 

RGS7 



MAF 

MYB 
NPRL3 
NFE2 
IL4 

ANPEP 

BACH2 

IRF4 

WHSC1 

NFATC1 



PITX3 

TH 

SLC6A3 

GDNF 

NR4A2 

MTA1 

FOXE3 

SLC18A2 

SOX2 

BDNF 

CHMP4B 



We then matched the candidate genes detected in PCG 
patients against their parents' CNVs to analyze the mode of 
inheritance. The identified CNVs were classified according to 
type (Figure 1) and we evaluated whether they were 
responsible for the disease phenotype of the patients or not. 

To narrow down the potential candidate CNVs (genes) 
and match the identified CNVs to target regions and/or genes, 
we first focused on known chromosomal loci for PCG, namely 
GLC3A(2p2-p21), which harbors CYP1B1, GLC3B (lp36.2- 
p36.1), and GLC3C (14q23). Second, candidate genes with 
identified CNVs were matched against both PCG-related 
genes and ASD causative genes including CYP1B1, MYOC, 
latent transforming growth factor-beta binding protein 2 
(LTBP2), PAX6, paired-like homeodomain transcription 
factor 2 (PITX2), PITX3, FOXC1, forkhead box E3 
(FOXE3), eyes absent 1 (EYA1), LIM homeobox transcription 
factor 1 beta (LMX1B), and v-maf avian musculoaponeurotic 
fibrosarcoma oncogene homolog (MAF) [11,16-18]. 

We next identified highly predicted functional partners 
of the target genes listed above using the Search Tool for 
Retrieval of Interacting Genes/Proteins database (STRING, 
version 8.3), and we determined if any CNVs were associated 
with these genes (Table 2). The STRING database comprises 
known and predicted protein interactions, including direct 
(physical) and indirect (functional) associations; these 
interactions are derived based on genomic context, high- 
throughput experiments, conserved co-expression, and 



previous knowledge. An example of the STRING 8.3 results 
for CYP1B1 is provided in Figure 2. 

RESULTS 

We performed aCGH to detect rare CNVs in PCG patients that 
were not present in 99 healthy controls. We identified 12 
patients with at least one rare gene-containing deletion or 
duplication, corresponding to 25 CNVs (5 deletions and 20 
duplications) at frequencies of 5%-30% in PCG patients and 
0% in the controls (Table 3). According to the literature, most 
CNVs have a Mendelian inheritance pattern [19,20]. 
Therefore, we performed aCGH using parent- offspring trios, 
because the use of family information can improve the 
sensitivity and specificity of CNV detection. We matched 25 
CNVs against the parents' CNVs, which were classified into 
three main categories: loss of 1 copy from one parent (5 
CNVs), gain of 1 copy from one parent (9 CNVs), and gain 
of 1 copy from de novo variation (6 CNVs). The remaining 5 
CNVs were classified as 'other'. No CNV was transmitted in 
an autosomal recessive inheritance manner. Considering that 
both parents of all patients were unaffected, loss or gain of 
one copy from only one parent is not likely to be associated 
with the disease phenotype. Only the gain of one copy due to 
de novo variation may represent an autosomal dominant mode 
of inheritance. 

We found that the 25 CNVs, identified within 28 genes, 
were not located in known chromosomal loci for PCG, namely 
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• CYP1B1 cytochrome P450, family l f Subfamily B, polypeptide 1; Cytochromes P450 are a group of heme thiol ate 



monooxygenasea. In liver mrcrosomes, this enzyme is involved i n an NADPH -dependent electron 
transport pathway, it oxidizes a variety of structurally unrelated compounds^ including steroids, fatty 
acids, and xenobioties (543 aa) 
(Homo sapiens) 

Predicted Functional Partners: 


"§ C * O 1/1 , 

a * s s e 

r c a u *j t; o . 

////// / iJ 


i ." COMT catechol- O-methyltra nsfera se : Catalyzes the O-methylation,. and thereby the inactivation, of cat [.,. J (271 a a) 
V GSTM1 glutathione S -transferase mu ij Conjugation of reduced glutathione to a wide number of exogenou [,.,] (213 aa) 

GSTPL glutathione s -transferase- pi L; Conjugation of reduced glutathione to a wide number of exogenou [..,] (210 aa) 
« UGTlAi UDP (jlucuronosyltransf erase l family., polypeptide Alj UDPGT is of major importance in the conju (...) (533 aa) 
'-" EPHXt epoxide hydrolase 1^ microsomal (xenobiotic); Biotransformation enzyme that catalyzes the hydro E--1 (*55 aa) 
^ GSTM3 glutathione s -transferase rnu 3 [brain)? conjugation of reduced glutathione to a wide number of (225 a a) 
3 UGT2B7 UDP glucuronosyltransf erase 2 family; polypeptide B7\ UDPGT is of major importance in the conju (..,) (529 aa) 
• UGTIAG UDP glucuronosyltransf erase 1 family, polypeptide A6j UDPGT is of major importance in the conju [--] (532 a a) 
■ft -GSTAL glutathione 5 -transferase alpha l ; Conjugation of reduced glutathione to a wide number of exege [...] (222 aa) 

UGT1A9 UDP glucuronosyltransf erase l family, polypeptide A9j UDPGT is of major importance in the conju £..,) (530 a a") 
















0.981 
0,977 
0.977 
0,977 
0.977 
0.976 
0,975 
0*974 
0.974 
0,973 



Figure 2. Evidence view of predicted functional partners of the CYP1B gene by STRING 8.3. Different line colors represent the types of 
evidence used to identify the associations. 



GLC3A (2p2-p21), GLC3B (Ip36.2-p36.1), or GLC3C 
(14q23). These CNVs did not include any target genes 
associated with PCG or ASD, nor the predicted functional 
partners of the target genes, and none of the genes had a 
specific gene function that appeared to be relevant to the 
pathogenesis of PCG (Table 3). 

These 25 identified CNVs might be PCG-specific 
pathogenic variants or may represent extremely rare benign 
variants that are not associated with the disease. Further 
genetic strategies to validate these CNVs are needed to 
identify specific gene functions relevant to the pathogenesis 
of PCG. 

DISCUSSION 

In this study, we performed aCGH in a series of 20 individuals 
diagnosed with PCG to discover novel copy number 
variations associated with this disease. 



First, we chose target genes by expanding the disease 
spectrum of PCG to ASD, resulting in inclusion of the 
following target genes: CYP1B1, MYOC, LTBP2, PAX6, 
PITX2, PITX3, FOXC1, FOXE3, EYA1, LMX1B, and MAP 
[11,16-18]. Unlike PCG, ASD is classified into different 
subtypes according to the features of malformation affecting 
the anterior segment structure, e.g., aniridia, Axenfeld-Rieger 
syndrome, iridogoniodysgenesis, Peters' anomaly, or 
posterior embryotoxon. However, some authors have claimed 
that PCG also involves abnormal development of Schlemm's 
canal and trabecular meshwork drainage structures. 
Surprisingly, mutations in the ASD genes sometimes cause 
PCG, and PCG genes can also cause ASD [2 1 ,22] . In addition, 
dysregulation or mutation of a few ocular genes can cause a 
range of clinical conditions. For example, PAX6 was first 
identified as the gene for aniridia, but is now known to 
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underlie a range of other ocular conditions including Peters' 
anomaly and a rare case of ASD [23-25]. PITX2 and 
FOXC1 mutations have been found in Peters' anomaly as well 
as PCG [26-28]. Peters' anomaly is also associated with 
mutations in two other genes, CYP1B1 and the FOXC1 -related 
gene, FOXE3 [29] . Therefore, we hypothesized that PCG may 
be considered part of the ASD spectrum; common genetic 
pathways may underlie these two disorders. 

We explored other genetic mechanisms underlying PCG 
by gene dosage analysis, because previous CYP1B1 gene 
mutation studies have been unable to explain the allelic 
heterogeneity and pathogenesis of PCG in all patients. 
Previously, Kim et al. [10] performed direct sequencing 
analysis of all coding exons and flanking intronic regions of 
CYP1B1 and MYOC in 85 Korean patients with PCG. These 
authors reported that about 70% of Korean PCG patients have 
neither CYP1B1 nor MYOC mutations (CYP1B1 mutation 
rate, 25.9%; MYOC mutation rate, 2.4%), results consistent 
with those reported for Japanese and Chinese patients. Indeed, 
12 out of 22 patients had only one mutant allele in the 
CYP1B1 gene [10]. In addition, the eye is known to be 
exquisitely sensitive to both reduced and increased gene 
dosage of key developmental genes. For example, gene 
dosage effects have been observed for PAX6 and FOXC1 in 
developmental ocular anomalies and ASD, respectively 
[12-15]. 

Overall, although we identified 25 CNVs (5 deletions and 
20 duplications) in 12 PCG patients, we were unable to 
correlate these CNVs with the pathogenesis of PCG using a 
reference-based approach. The identified CNVs were not 
located in known chromosomal loci for PCG, namely GLC3 A 
harboring CYP1B1 (2p21), GLC3B (Ip36.2-p36.1), or 
GLC3C (14q23), and did not include any target genes 
associated with PCG and ASD, nor highly predicted 
functional partners of target genes. However, our data suggest 
that altered gene dosage may explain the disease pathogenesis 
in PCG if these CNVs are PCG-specific pathogenic variants. 

Our study has some limitations. We were unable to 
exclude the existence of pathogenic variants that were too 
small to be detected using our platform (<500 bp). Despite the 
high genomic resolution of aCGH used in our study, some 
genomic regions might not be covered well or may not have 
been assessed due to technical difficulties. Furthermore, we 
did not perform further studies to determine the gene functions 
or expression levels of the CNVs we identified. 

In conclusion, this is the first study to comprehensively 
investigate gene dosage effects in PCG. We believe that the 
preliminary results and the CNV filtering strategy that we used 
can broaden our understanding of the genetic mechanisms 
underlying PCG. 
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